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EIGURES 



iSSTMOT 


This study concerns with, the inviscid hypersonic 
flow over plane and axisymmetric hlunt bodies producing a 
power-law shock. Unifornil 3 '' valid solutions, far downstream 
from the blunt nose, have been obtained for the inverse 
problem, vi/hich prescribes the shock leaving the body to 
be determined, using the method of ’inner and outer 
expansions' originated by Van Dyke. 

The significance of the entropy layer for various 
power-law shock has been discussed. Power-law bodies for 
various values of the exponent, n have been constructed 
for the axisymmetric cas The second-order results for 
n = 1/2 and 2/3 have been obtained. First-order body and 
flow variables ha.ve been obtained for n = 2/5 in the planar 


case 
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u 


constant, (3.1) and ( 4 .I) 

constant , 2 

constant, 

density in the outer region 

similarity variable for velocity, v 

similarity’- varia-ble for density, 

similarity variable for pressure, p 

Integer specif y’-ing number of dimensions 
= 0 planar case, 

= 1 axisymmetric case 

pov/er-law exponent 

(l-n)/n 
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2(l-n) 

pressure in the outer region 
pressure in the inner region 

velocity vector 

density in the outer region 

time 

velocity in the outer region parallel to the 
x-direction 

velocity in the inner region parallel to the 
x-direction 


V 


V 
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A 


cr 

M 


velocity in tlie outer region parallel to the 
y~direction 

velocity in the inner region parallel to the 
y-direction 

coordinates along the body axis 
coordinates in the lateral direction 

proportional distance across the shock layer 

density in the inner region 

ratio of specific heats 

stream function 

local shock angle 
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Subscripts : 

b evaluated at body surface 

1 first order 

2 second order 


CO 


free stream conditions 
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GHiil’TER 1 


nT JRODUGTIOK 

Hypersonic flow past iDlimt bodies v>?ith sboclc waves - of 
very gree^t intensity tliat h.a.vo a, power-law form are considered 
for the plane and smi-sjaametric cases, for the power-law 
exponent, n < 1 , the shock shape being blunt near the nose, 
the stream lines pa,ssing through it a,re subjected to large 
changes in entropy resulting in the formation of a, layer of 
rela.tively high entropy g-as. This layer sometimes referred to 
as the 'Entropy Layer' exhibits flow properties quite similar 
to a boundary layer and is different from the slightly disturbed 

I 

outer flow field'. The influence of this entropy layer on the 

3 7 

asymptotic flow field horn been analysed by Sychev, Yakura and - 
several others. 'The outer region where the streamlines cross 
the oblique portion of the shock-wave can be treated by hyper- 
sonic slender body theory ( .1 ) , Sychev has considered the 
inverse problem of determining the body that produces a power- 
1, aw shock-wave specifying a paraboloidal shock-wave and properly 

accounting for the entropy layer he has shown a method to con- 

7 

struct the body shape nimerically, Yakura, has considered the 
inverse problem using the method of inner and outer expansion 
(8 ) to obtain uniformly valid solutions far downstream from 
the nose of a sli^tly blunted body in hypersonic flow. An 
attempt has been made in this work to construct aiiaJ-ytical 
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solutions for the inverse prohlem, prescribing a power-law 
shock-wave. The significance of the entropy layer for various 
powcr-lavr exponents have been discussed and power- law bodies 
constructed. 

For the slender bodies the line of development has been 
concerned with the small disturbance theory, where it is 
assumed that the velocity perturbations are small compared with 
the free stream velocity. The velocity perturba.tions, however, 
are not small compared with the loco.1 sound velocity and the 
pressure perturbations are not small compared with the free 
stream static pressure. Hence the disturbances are not at all 
small in the sense usually associated with the linearised super- 
sonic theory. This theory started with the hypersonic similitude 
principle. The idea is due to Tsien (1 ), thou^ Hayes introdu- 
ced the general hypersonic similitude principle which states 
that a small disturba.nce hypersonic flow' is equivalent t.o an 
unsteady flow in one less spare dimension. Particular group of 
solutions to the small disturbance equ— ations that have been 
v/orked out have been almost exclusively of the self-similar 
typo. They arc generally termed similar solutions, and have 
the property that the solution in the lateral variable y or r 
at one value of t (or of x for the corresponding hypersonic 
steady flow) is similar to the solution at ary other value of 
t. This property permits a decrease in the number of indepen- 
dent variables from two to one thereby reducing the partial 
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differential equation to the ordinary one, the treatment of 
which is much easier. Many of the investigations of self simi- 
lar solutions of the blast wave type 8,re directly applicable 
to the hypersonic steady flow problems, Lees and Kubota have 
used the blast wave theory for simula,ting the flow past blimt 
bodies. But the a,nalogy between the flow field of an intense 
explosion and the steady transverse flow field of a blunt body 
moving hypersonically in one more space dimension is a small 
disturbance approximation va.lid only fo,r from the body, The 
strong perturbations due to the blunted nose creates a region 
of high ent 2 ?opy gas near the body in which the small distur- 
bance equations and the hypersonic similitude principle fail 
to apply, Iho entropy layer is similar to the boundary layer, 
in that the pressure is essentially constant across it and the 
velocity gradients are large, 'The entropy layer also has a dis- 
placement effect on the flov/, like the boundary layer, 

3 

Sychev considers the inverse problem of determining the 
body that produces the outer small-disturbance flow field given 
by the analogous blast- wave. He recognizes the shortcomings 
of the blo.st v;ave theory, and for the case of plane and axi- 
symmetric gas flow v^ith shock waves of power- law form, he shows 
that the use in the flow problem of the exact solution for the 
corresponding unsteady self similar gas motion requires a 
supplementary refinement of the thickness of the high entropy 
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layer. He also sliows a methocl for introducing such, a correc- 
tion and constructing the shape of the body contour on which 
is to he applied the pressure distribution on the ba,sis of 
the theory of small disturbances* From his study he concludes 
thc,t ; i) the 'equivalence principle’ or the law of plane 
sections is not valid in the entropy layer, ii) The blast wave 
analogy applies t<'' a specific body that has a relatively large 
\ lateral dimensions compared with the nose dimensions. nO In 
the range of values of,*/ of practical interest there exists 
a definite interval of values of the exponent n, 2/3+d < ^ 

< n (o = 0 or 1 for planar or axi-symmetric cases and 
n* = (2 + 2/’^ ) / (3 + #+ 2/Tf?)) where proper consideration of 
the entropy layer is necessa.ry ■ for e.g. in determining the body 
contour for a given shape of the shock-wave. Specifying a 
- parabolic shock-wave and properly accounting for the entropy 
layer, he calculates a body shape that becomes very large far 
downstream. And finally he concludes that the self similar 
motior^in all cases certainly retain their significance as 

asymptotic representat-tons of the exact solution. 

u. 

It was first shown by Guirc^ in his study of the flow 

past a blunt flat plate that the entropy layer is a region 

of non-unifoimity , characteristic of a singular perturbation 

problem and can therefore be treated by a perturbation method 

Q 

first applied to the bounda,iy layers. Guirad concludes that 
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the blast wave theoiy is correct as far as the leading term 
is concerned, the entropy layer being felt only through the 
first correction, thereby settling the point which has been 
controversial in the past namely, whether the displacement 
effect of the entropy layer generated by the detached shock, 
ahead of the nose, Invalidates the blast wave analogy theory. 

Yakura uses the method of inner and outer expansion in 

obtaining uniformly va,lid solutions far downstream from the 

nose of a sli^tly blunted body in hypersonic flow. This method, 

8 

originated by M.D. Van lyke is applied to the inverse problem. 

He transforms the equation of motion such that the outer region ^ 
of the shock layer a,nd the inner region corresponding to the 
entropy layer can be treated separately. With a proper choice 
of variables he forms the asjonptotic expansions for each region.; 
Alth 'Ugh real gas effects become significant at hypersoiip speeds: 
he restricts the :enalysis to the case of an inviscid, non 
heat conducting perfect gas, and to the limiting cane of Mach 
number tending to infinity. He treats the plane and axisymme- : 
trie flows past bodies producing hyperbolic and power-law shocks; 
The iiyperbolic shocks correspond to blunted wedges and cones ' 
in two and three dimensions respectively. The exponents of the i 
power-law shocks considered are n = 1/2 for axi— symmetric case 
and n = 2/3 for the plane case. Yakura concludes that the [ 

entropy layer is, in many respects, analogous to Prandtl’s 
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viscous boundary la,yer with, concion features of displacement 
thickness, zero norma-1 pressure gra,dient, and a velocity 
gradient. The flow field influenced by the entropy layer varies 
with body or shock geometry. When n = 1 (wedge or cone) the 
effects of bluntness become (for large x) concentrated near 
the surface; whereas for n = 1/2 they spread out into the flow 
field having the same rate of grovrth as the body. Thus, the 
relative thiclmess of the entropy layer (compared with the 
body radius) becomes appreciable only for values of n much 
less than unity. 

Presented in this study are anal 3 rtical solutions of 
the flow field in the form of asymptotic expansions uniformly 
valid for largo distances downstream. The examples treated in 
this study are plane and axi-symmetric flows past power-law 
shocks of any exponent n(2/3+3 < n < 1). Yakura^s method is 
extended to the general case of power-law shocks. Here again 
the equations of motion are transform.ed such that the outer 
region of the shock layer and inner entropy layer can be 
treated separately. With a proper choice of variables, asympt- 
otic expansions are formed for eo.ch region. The outer limit 
process transfers the region of non-uniformity (entropy layer) 
to the inner unknown boundary (body) and the outer expansion 
formed is valid every where except at the inner boundary. The 
inner limit process transfers the outer region to infinity 
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and the in., er expansion is valid in the inner region hut not 
far away. The limiting process is equivalent to shrinking the 
co-ordinate system or physically moving to infinity in the 
dov/nstrean direction with the nose hluntness fixed. It is the 
inner expansion which deteimines the body. Erom a physical 
view point, the outer limit consists in moving far downstream 
(or, alternatively lotting the nose radius shrink to zero) 
while remaining away from the surface, whereas in the corres- 
ponding inner limit one stays v.’ithin some multiple of the nose 
radius of the surface. Boundary conditions for the outer problem 
are the shock .conditions and for the inner problem use can be 
made of matching between the inner and outer expansions to 
extract the inner boundary conditions. A uniformly valid solu- 
tion is then constructed by forming a composite expansion 
consisting of both the inner and outer exp.ansions . 


1 ! 



CHAPTER II 
GEHERAL^ 

Plow Equations : 

Steady flow of an inviscid perfect gas with constant 
specific heats is described by the following set of equations, 


div (^q) = 0 ■ 

(2.1a) 

^ (q « grad)q'+ grad p = 0 

(2.1b) 

q-grad (p/pV) = 0 

(2.1c) 

The conservation equations written 

cartesian co-ordinates are 

in 2-diinensional 

(pu) + ly ((Jt) = 0 

(2.2a) 

r („ 0Uv ^ _ P 

V 9x ^ ^ ^ 0X “ ^ 

(2.2b) 

9v , 9 vn , 0p p, 

^ 9y = ° 

(2.2c) 

u (p/pt) + ■^ ly (p/f^) = 0 

(2. 2d) 


Equation (2, 2d) expresses the fact that the entropy 


is constant along streamlines between the shockwaves and 
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the body. I’his property has been exploited by defining a stream 
function that satisfies the continuity equation as an indepen- 
variable. usine Yon-Mises transformation. Then the differential 


equations (2.2a - e) become 

^ = Z 

dx u 

2 2 2-f , 

u + V + p/p = 

p/^.-y = f(l|r) 

where j = 0 for plane flow and 


(2.3a) 

(2.3b) 

■(2.3c) 

1 (2,3d) 

(2.3e) 

i = 1 for axi-S3nmnetric flow. 


Bound ary Oo nditio ns % 

Blow variables have been made dimensionless by 
referring velocities and density to their free streajm values 
u^ and and pressure to . Equations (2.3a-e) are 

unchanged with respect to this transformation. Boundary 
conditions are the flov\/ quantities behind the shock-wave and 
are obtained from the oblique shock relations and expressed 
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in tlie form for M = °° 


^ = 1 + 7 sinV (1 - e) 

Jr CO 


where 


'<-1 

V+1 


oo 


2 

"Uoo 


Pc 

^1" 


Dividing throughout by '/ Iv![^ v/e get 


r u 

oo oo 


1 2-2 
sin^cr 

7 X ^ 


+ 1 


'Ko "*■ °° we get the dimensionless pressure as 


2 . 2 
p = ^rT~T 


p "^+1 

— = or in the dimensionless form 

P :c ^ 


p - 1±™L 

u = cos O’' 

^ / 

V = 6 sin cr 


where u and v^'are the velocities tangential and normal to the 
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shock and o~ is the local shock angle. 


2 ~f— 1 2 

u = cos cr + sin cr 

1 +1 


2 . 2 
1 ~ ~ sin cr 


Y = COS cr sin cr 


y-1 

-f+T 


cos crsin o" 


Y +1 


cos cr sin cr 


y = [ (1+d) f ] 


V(l+3) 


Boundary conditions at the shock are 


P = TTY sin cr 


1 ,± 1 . 

Y-l 


u 


2 ...2 

-T+T 


1 - vr-T- Sin O' 


V = sin cr cos <T 


J = [ {i+ j ) f ] 


V(1+3) 


(2.4-a) 

(2.4h) 

(2.4c) 

(2.4d) 

(2.4e) 
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Shock-Wf'V 6; 

Tlae blunt shock-wave is taken to have an analytic 
functional form 


y 


Ax 


n 


y is the normal co-ordinate with respect to the free 
stream direction x, in both plane and axi-syimnetric flows. 

The local shock angle is given by 

n— 1 

tan cr = An x 

The form of the stream function at the shock is given 
as 

Ijr = (2.5 

1+D 

Equations (2,3a-e) are the set of five equations in 
the five unknowns, (p,o, u, v and y) as functions of the two 
independent variables x and l|J subject the shock conditions 
(2.4a-e). Equation (2.3©) gives an insight into the extent 
of T:hc entropy layer. The entropy function p/^’-'" on any stream- 
line depends only upon the local shock angle evaluated at the 
point of intersection of the shock and the stream-line. The 
sti-eamlines which intersect the shock far downstream and far 
out from the body define an outer region wfhere the stream- 
line is of the order of In this region the entropy 
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function takes on its small disturbance value or approaches 
a Viniform value depending on whether n < 1 or n = 1. The region 
in v/hich ijJ is of the ordei’ of 1 corresponds to the entropy 
layer. In this ^inner region’ the stream-lines origina,te in 
the vicinity of the nose and p/p"^ is relatively large. The shock 
is assumed to bo smooth and so the entropy is a continuous 
function of the stream function. Tor thorough investigation 
of the two regions of interest the method of inner and outer 
expansions have been used. 

I nner and Out er Expans ion s 

Asymptotic expansion for large down-stream distances 
have been developed. The small perturbation parameter associated 
with the limits and expansion is 1/x where x is the longitudi- 
nal down-stream co-ordinate as shown in Fig, 1. A new indepen- 
dent varial;le of ord.er unity is introduced in the outer-region 
extending to the shock. 

Outer Independent Variables 

X = X 

w = 

Outer Dependent Variables 


P, D, Uj V, Y which are respectively the 



14 


pressure, density, velocities along and perpendicular to tlie 
X-direction and the normal co-ordinate, 

'This choice of the outer independent variable is fully 
consistent with the definition of the outer region, where if” 
is of the order of 

Let P(X,w) denote pressure in the outer region or in 
general let pCXjw) he any flow quantity in the outer region. 
Outer limit is defined as the limit as the perturbation para- 
meter tends to zero with the outer variables fixed i.e, 

0 lim []P(X,w)]= lim FCXjw) 

as X keeping w fixed. 

In the outer region w is of the order one. In the 
inner region f is of the order one and the inner independent 
variables are x and ijl and the dependent variables are p, ^ s 
u, V and y. The inner limit is defined as the limit of a 
flov/ quantity as x keeping T}r fixed, let f (x,l|r) denote 
a floo»/ qpaantity 


I lim [ f (x/[jr) j = lim f (x, f) 

as X -* keeping ijJ fixed. 

The outer limit consists in moving away from the 
surface some appreciable fraction of the distance to the shock 
and the inner limit considers the region near the surface. 
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The eauations of Section 2,3 are transformed into the new 
iriiier an;; outer variables. The ’outer equations’ y\/hich are 
functions of (X,?;) are valid in the outer region and they will 
be used to determine the outer small disturbance solutions, 

Tor the outer equations transformed shock conditions are the 
boundary conditions. The inner equations are dependent on 
(Xjljr) and are restricted to the inner region and these deter- 
mine the inner solution taking care of the displacement effect 
of the entrop 3 ^ layer. Inner equationsas a result of the asym- 
^ ptotic expansion suffer a loss of boundary conditions. 

Composite Ex pansion 

Composite expansion is constructed from the inner 
and outer expansion and is valid in the entire flow field. 

Both the inner and outer expansions are not valid in the 
regions of their counteirpart , The sum of the inner and outer 
expansions aro taken, outer expansion being rewritten in the 
inner variables and the enter expansion of the inner expansion 
is subtracted to remove duplication. The composite expansion 
takes the form 

f (x,ljF) = I €x[ fCx,'^') ] + 0 ex [ T(x,ljr) ] - 0 ex 

[ I ex[ f{x,f) ] ] 

Near the body, as w 0, contributions of the outer 
expansions cancel out. 
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;^o ck--Wav8 of the P ow er- Law Pqpn 

AS already stated the theory of small disturbances 

in a hypersonic stream, the problem of flow past a plane or 

accisymmctric shock-wave or body is equivalent to the problem 

of one-dimensiona,l unsteady gas motion under the action of a 

plane or cylindrical! piston. In this analogy the class of self 

similar motions with very intense shock-waves propagating 

according to a power-law corresponds to a class of steady 

n 

flows with shock->vaves of a power law form j = Ax. . With the 
mach number of the undisturbed stream °° values of the 

exponent lying in the interval 2/3+ j < n < 1 , where 3 = 0 or 1 , 
correspond to flows past convex bodies of pov\/er-lav/ form. 

The case n = 2/3+d is singular and corresponds to the problem 
of a strong, explosion. In this case the interpretation as a 
flow problem consists in the assumption of a finite drag force 
acting on the leading edg;e of a body of vanishing thickness. 

In other words there is an analogy between the appearance of 
a strong explosion and the effect of blunting the leading edge 
of a slender body at large distances from the bluntness. 

At sufficiently large distances from the leading edge 
of the body the transverse velocity component v is propor- 
tional to the local angle of inclination of the shock surface 
*0" ' and the pressure to the square of the sine of the angle. 


■r 
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^ . P ~ 0-2 , cr ~ 

X 

Substituting these estimates in equation (2 ,3c) and (2,3^) we 
find that the relative change in pressure across the entropy 
layer is 


AE ^ (T ^(1+d)/(1-ii) 

P 

. < 0-2 for n > 2/3-f j 

5 

Since the small disturbance theory involves just the 
same limit of accuracy the change in pressure can be neglected, 
i'or the estimate of the relative thickness of the entropy 
la3’er Equations (2.3e) and (2,3b) are used. Erom the condition 
of constancy of the entropy along stream-lines and the boundary 
conditions for p and p we find that along the entire entropy 
layer 



Using the estimate for the pressure p, we obtain the following 
estimate for p 

Taking u<'-U^ using (2,3b) we find that relative thickness of 
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the entropy layer is 





n(1+j)/(l-n ) 


- 2/4 


It is negligihly small whe^ithe exponent of cr is greater than 
2 i.e, 

n(1+j )/(l.~n) - If >. 2 
n > n* = (2+2/'/)/(3+o+2/-/ ) 

This means that in the range of values of "( of practical 

interest there exists a definite interval of values of the 

* 

exponent n, 2/3+3 £ where proper' consideration of the 

entrop 3 ’' layer is necessary for e.g, in determining the body 
contour for a given shape of shoclc-wave. An attempt has been 
made in this present work to find out the effect of entropy 
layer for various n. 



CHikPTER III 


BLm^r AZISYilviE'IRIC CASS 


Tnc two different cases name!;/ the axisjcmetric and 
plane flov/s pive rise to different "boundary conditions and 
hence different outer expansion. Therefore the tv\/o cases are 
treated separately. 

The shock in Sig. (1 ) is given "by 

Y = A '(3.1) 

Prom equation (2,5) the form of the stream function is 

f = [ ]^/2 ( 3 . 2 ) 

The basic equations for axi symmetric case are 


^ ^ Z 

dx u 


(3.3a) 


y 


if = 



(3.3b) 



(3.3c) 


2 2 2Y p 

u + V + ^ 


1 


(3.3d) 




(i=l) .1^ 

''Y+l ^ '/+1 


. 2 
srn cr 


(3.3e) 
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sin cr and cos o" are found in terms of stream function f” as 



sin^o" = 

(3.4a) 





n*. 

2 Bt 

COS (T = Jr_ — 

1+Bijr^'' 

(3.4b) 

where 

1— n j 

n-| — and 

O/f ) 2/n 2 

B = 2 / A / n 


Boundary conditions follow from (2.4) 



p = 2X-^+D (1+BllF ') 

(3.5a) 


p = ('^+1)/('f-l) 

(3.5b) 


u = 1 - 2/(-/+i)(i+Bijr ) 

(3.5c) 


V = 2 ■, ifBf ^ 

(3.5ci) 


y =-j2’i 

(3.5e) 

Outer 

Exp alls ions 



The outer variables are chosen as 

(X,w) where 


X = X 

(3.6a) 


II 

(3.6b) 
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The outer expansion is valid in the region avi/ay from 
the body and w is so chosen to be of the order one so that the 
outer variables are expanded in pov;ers of x for large x 
distances . 

Equations (3.3a-e) are transformed into the following 
equations ’.vith P, D, U, V and Y as the dependent variables 


Y ^ - IF 4. 9-n 
^ ax “ u ^ Dtr 


T il - y 
^ aw “ DU 


2n 


2n av . 2n-1 ^ 

^ ^ - 2n X ^ ^ + Y ^ 


0 


„2 ,.2 2Y P . 


P . _ .//-I y 

“ V+1 ^ Y+1 


1 

1 +B (wX^^)^1 


(3.7a) 

(3,7b) 

(3.7c) 

(3.7d) 

C3.7e) 


The shod conditions are transformed into the following 

P >2_ 1__ 

■f +1 n^ 

1+B (wX^^) 

= 1 r 1 - + 

y+1 n^ no L n. n 

w ' BX Bw ' X ^ 


] (3.8a) 
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l! 


T+i 

■y-1 


(3.81d) 


U = 1 - 


1 


•f'+T "“n n 
w 'BX ^ 


[ 1 — h-^ + ] (3-8c) 


Bw ‘X 


V = 


2 1 
y+1 — 

w ^ X 2 


[1 - 


1 


^1 ^2 
Bvv ^X 


+ ...](3.8d) 


I Si 

Y = X 2 wX 


. ( 3 . 8 e) 


at 


w = A /2 


v/liere 112 = 2n 


' 2 lie asymptotic expansions for the outer flow quantities are 
Ghe following;, v/hich can be arrived at from equations (3.8a-e) 

^-1 (w) 


l-(X,w) = 


n. 


+ 


2n. 


+ • . • 


X 


X 


{3.9a) 




Do (w ) 

•D(X,w) = I)^{w) + -4: + » 


n. 


X 


(3.9h) 


U^w) U^Cw) 


IT(X,W) = 1 + -5- + 

X X 


H” * * 


(3.9c) 


4(w) Y^(v/) 

Y(X,w) = + “3S72 + ••• 

X ^ X ^ 


(3.9d) 
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f(X,w) = Y^(w) X + + ... 

X ^ ‘ 


C3.9e) 


Substitu'l^ing the above expansions into equations (3,7a--e) we 
get the following equations. 

First Outer Equation: 


Y 


2w 


1 ^ 


Y. 


n 


(3.10a) 


. 1 
■1 dw - D, 


(3.10b) 


V + ^ 

^ 1 n^ dw 


2 '^^l 

5^ ^1 ai“ = 0 


(3.10c) 


D 


2 

W+r q + 1 


■L 

n. 


Bw 


(3.10d) 


+ Y, + 


Tr2 , 2'/ ^1 


1 ^ — “TXj" = 0 


(3.10e) 


Boundary conditions 
P 


1 “ V+1 


2 aV 


(3.11a) 


D - -i±i 

■^1 “ 7-1 


(3.11b) 


2 2 2 

u. ^ A n 

^ y+1 


(3.11c) 
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(5.11d) 


Y. = A 


(3.11e) 


Second Outer Equations; 


. U. ■ D- 

nY^Y^ + (n-n2) Y^Y^ - + J ^ + Y^Y^ + 2rav [ - ^ _ >2^ ] 

1 


(3.12a) 


dYp dY. 

^F" + l2d?r 




(3.12b) 


Y + "^2 2 r „ 

^2 + 5n7 dF^ ~ — — [ Y 
1 3no 


dPp dP. 

1 3r-+ ^2 35r ] = 0 (3.120) 


2 P 7^ ^ 9 ^ i 

2^2 ^1 ■*■ ^^1^2 FF t PT ^ — ] = 0 (3.12e) 


1 D. 


^1^2 


^ Y+l ^ Y+1 


B w ■ 1 


(3.1 2d) 


Boundary conditions 


P - > . 2__ 4 4 
■^2 “ Y+1 


(3.13a) 


0 


(3.13b) 
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(3.15c) 
(3.13d) 
(3.13e) 

Solution of the Outer Equation: 

Equations (3.10a,'b,c and d) are to be solved for the 
unknowns Pj P, Y, and Y and U can be calculated from equa-tion 
(3.10e). Two different methods have been adopted for the solu- 
tion of the first outer equation^ namely^ the similarity approach 
and the numerical integration by the Runge-Kutta method. 


Similarity Approach s 

The results of integration of 

similar gas motion (cylindrical case) 

the relations 

the equation of self 

appear in the form of 

P = Pg(x) h(A) 

(3.14a) 


(3.14b) 

V = v^(x) f(^) 

0 

(3.14c) 

X = y/y^ 

(3.14d) 


where the index s refers to the condition on the shock-v/ave. 
The details of this method are given in Appendix I, The 


^^2 - + “TkT 




A^n 


TT ■ _ 2 .3^3 

^2 - - A n 


Y 2 = 0 
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functions f (/*), g(^-) and ii(A) are a.lso described in Appendix I, 
‘The tv/o term outer solution hecomes 


P(X,w$ = Pg(X) la(k) 


2 

i+f 


[ A^-n 


2„2 1 


X^2 


R(X,w) = Rgg(A) 


A —I 

-^ + 1 


x2^2 


] 


(3.15a) 


= g(^) 


(3.15b) 


Y(X,w) = V^(X) f(>0 


= [ 


2 1 
iV\' ^ n /2 
X ^ 


+ 1 




X 


(3.15c) 


Y(X,w) = A AX' 


n 


(3.15d) 


Relation between w and follows from equation (5* 10a.) 

,2 


V'/ 


= [ ('(+'* ) ^ f (^) 1 


(3.15e) 


U velocity is determined from equations (5.10e) and (3.1 2e). 
Two term outer solution for the velocity turns out to be the 
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following: 


U(X,w) 



1 


" ^ 


[ (- 


2 y 


aV( 




2 


+ 


8 


(-/+1)' 


/. 4- ^^2/v> 
A n f (A; 


+ 


2y 


(- 


2(i^~1) 

(y+1)^ 


A‘^n'^li(>') ) ] 


1 


2no 
X 2 


(3.15f) 


numerical Solution (for the first Outer Equation); 

The following differential equations are formed from 
equations (d.lOa-e) for the dependent variables 
and . 

dP. 

35 ^ = Pf = Pf/Pfg (3.16a) 


Pf^ = (nY^^ - 


w 2w v^/Y 


1 


and 



A*)/-^ “VT- 


where 



dD^ 

dw 


(3.16b) 


Pfw ^ 


/ w 


dV^ 

dvT 


2Y,Pf n. 


dY, 

"W 


DiYi 


U, = 


V 


X.!i 

V-1 


(3 . 1 6c ) 


(3.16d) 


(3.1 6e) 


Range-Kutta integration is carried out using the 
conditions at the shock from equations (3.11a-e). 


Inner Expansions; 

I" and X are the independent variables and equations 
(3.3a~e) describe the inner flow field, Ihe outer solutions 
are given as functions of A from the self similar solution. 

So the inner expansion can-^not be arrived at by studying the 
behaviour of the outer solution as w — 0, Hov/ever, if we 
investigate the asymptotic behaviour of the outer equations as 
w becomes small, the inner expansions can be arrived at. 

Prom 3.10c and 3,12 c we find the asymptotic behavioiir 


of P(X,w) to be 
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P(X,w) 


P^(w) P 2 (w) 

X^2 x^2 


(3.17a) 


Where 


1 


h ^2 

(w ) = ^ w 


+ p 


It) 


P2(w) 




3 112 

— w + P2T, 


where the subscript ( )-j^ refers to the flow wariables 
evaluated at the body. So the inner expansion for pressure 
can be written as follows in terms of inner variables. 


p(x,f) 




(3.17b) 


Prom similar considerations the expansions for the other 
variables can be written as follows; 


Prom equa.tions 3.10d and 3.1 2d we get 

^o(w) 

D(X,w) = D^(w) + 


(3.17c) 


D^(w) = [P^(w)] ^ (^) 
I’o (w ) 

^ 2 <^> = 


where 
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From 3.10a and 3.10e we get 


Y (w ) — — — ^ rr-r— — i— 

(3.17d) 

(w) = nY^(w) 

(3.17e) 

- TTr-s:^ 

(3.17f) 

Similarly from 3.12a and 3.12e we get 


^ V. (w)U.(w) 

Y^Cw) = 

(3.17g) 

V^Cw) = nY2(v^) ' 

(3.17h) 

ord 

U2(w) = V^(w)V2(w) 

(3.171c) 

From rlie above relations the inner expansion is found to be 

the following: 

/ ^ PiCf) PpCt) - . 

p(x,f) _ + — 2nr> *** (3*18a) 

x^2 x^2 

flCf) / V 
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u(x,i|r) 


v(x,ilr) 


y (x,f) 


u^Cijr) 

2 (n+n^-n^/Y ) 


-y-l (1l) 

~ 2f-n2 ( 1 -n^/ 2y+ii^-n^/ ) 

“~2T'" X ^ 


nV2/ 3np/2'/ - 

(f )x + y2(f)x; ^ + ... 


(3.18c) 


(3.18d) 


(3.18e) 


Substituting tbese expansions in tbe equations (3.3a-e) we 
get the following set of first and second inner equations. 


First inner equations: 


n<^ -[ 

■“ 2 “ 


ir. 


(3.19a) 


^1 


d^ 



(3.19b) 


dp. 

dr= 0 


(3. 19c) 


ui+ 






h 

?i 


= 0 


(3.19d) 


Pi .Y-.1^^ 2 1 

T^T" (l+Bfi^i ) 


(3.19e) 
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Second inner equations: 

Up 

^2^2T “ ^2^ ^ ^2 (3.20a) 

dy. dyp 

^Sd^+^ldf (3.20b) 

. (3.20c) 

1 / 2 ) , 

(3.20d) 
(3.20e) 

lirst Imicr Solution: 

It is clear from the first inner equations (3.18a-e) 
that tv;o matching conditions in p and y are required from the 
outer solution. The following are the conditions written in 
inner variables. 

1 

X^2 


U 2 = 0 


Pg = 




dpp 

= 0 (when n / 1/2) 

dpp V. ^2“^^ 

dlf ' y^ ^ 2^7 ^ (when n = 


P-| (^r1^) — ^ n^( /-|^) 


(3.21a) 
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(3.21b) 


V;lien tiie outer solution is numerics.lly integrated tbe values 
of and at w = 0 are taken .for matcblng with tbe inner 
expansions* 


Then the first inner solution is the following: 
2 .2 2 


Pl(x,ir) = A n 

X 2 


(3.22a) 


f^(x,ltl) = |^(AVh(A^)) (1+Bf ^) — 4-^77 


x""2/ 


(3.22b) 


U.| (x,l(j) = 1 


2Y (A^n\(A^)) 


f-1// 


/ X 




(3.22c) 


v^ (x, |J) - ^|j(2.il) [ -72-4 


A n h( A-j^) 


] 


I/ 2 A 1/2 



p-nyif.' 


+ 


X 


^b 1 

T=5:J^2V| 


(3.22d) 


y.,(x4) = 


'/ i V2 

(2.^P(f)) (■ 


aAA ( x .|^) 


1/2'/np/2/ n 72 -/ 
“) X ^ + A A, X 


(3.22e) 
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v/here Pdfj) = / 


1 

(1+Bljj^l)"*^'^ 



Composite Expansion: 


ihe first order solution is given by the construction 
of the composite expansion written in inner variables as 
described in Chapter II. 


p^(x,i|r) = 


x^2 


(3.23a) 


xV '^ 

(3.23b) 




C'+1)‘ 


(1+Blf^l)''/^ (BlfPl ) 


2y 


(1^+1)‘ 


A 


1 -.1 


x^2 


(3.23c) 




2Aa f(A) 


n. 


/-nV2, 


1/2 


1/2 




A 


■n^( A, )' x^^^^2^ 


(3.24d) 
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Tl(x,ljr) = , ,1/2/ 




1/2 

((P(i}r)) -p^(ijr))] 


(3,23e) 


where 


(1J) 


= ht / 

]|j oo 


dt 


Second Inner Solution: 

Here again the matching conditions are 

PjCx.lT) = - ^ A'^n'^hC X^) 


(3.24a) 


= 0 


(3.24h) 


From equation (3.20h) 


^2 = T7 


(3.25a) 


Fp(llr) = - 


/ 2 2 . ^ 

av A^n^( X^) ^ 

(A^n^(A^) ^'^(1+Bf^'')^^ 


P 2 = ~ fipf A^n'^hC A^) (when n/1/2) 


(3.25h) 


v/here 
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n. 


(2'P 


(n2“2-/ 


2 , 4 ^ 4 ,^ 




(3.251).) 


(when n = 1/2) 


Yr 




f 


2 


■2 v'i 
■V P-i 


(3.25c) 

(3.25(3) 


Ug = 0 


(3.25e) 


She second inner solution can be written directly from 
(3.25a-e) 


2np 

P2(x,l|r) = Pg/x 

(3.26a) 

n./''' +n. 

4p(x,l|J) = p^/x ^ 

(3.26b) 

U2(x,ljl) = 0 

(3.26c) 

(2Y-n2)/2f + {rirf-nJi ) 

Y 2 (x j ^ ) = Y 2 / X 

(3.26d} 

y2(..t) = 72 ■ "2 

(3.26e) 


Ihe composite expansion for the second term can be 
found out in the same manner as before 
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where 


P 2 (x,l|J) = - k^n\{A)/x^Z (when n ^1/2) (3.27a^) 


2n 

PgCxjljr) = - [ A^n\(A) + (n 2 - 2 '/' )l{r ]/x ^ 


Cl 


(when n = 1/2) (3*27v') 


(x,|r) = p (xjljJ) 






(3.27b) 

2 


(V+1) 


+ tI-T AA4f2(A) + 


(- ^A^n^h (X) l/x ^ 

(■/+1 r 


(3.27c) 


= - “pj- A^n^ -g- ' 72 f (X) + V2(x,f) •■ 

x 2 


n 


2 (27-n^2Y+ (n^-nV/) 

(^ - npFjWA ^ ^ ^ 


(3.27d) 


, ,, , (2V-np)/2/ + (n^-np//) 

y2(x,ljr) = ~ 5'3(f)A ^ ^ 


(3.27e) 


F, (1r) = Lt ^ / Pg 

It oo 1 
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The complete second order solution follows from equations 


(3*23a--c) and (3.27a~e) 

p(z,l|l) = p^(x,lll) + (3.28a) 

= p^(x,Tjr) + (3.28b: 

u(x,f) = u^(x,^) + U 2 (x,l|r) (3.28c) 

^(xjljr) = w^(x,Tjr) + V2(x,f) (3.28d) 

y(x,l|r) = yj(x,l|r) + y 2 (x,l{r) (3.28e) 


The results of the second order solution are given 
in figures ( 8 -ij ) , 



CHAPTER 17 


P LAICE ppy/ER LiW shoce: 


The shock is given by 

y = Ax^ (4.1) 

Prom Equation (2.5) the form of the stream function is 
given by 


I|r = Ax^ 


(4.2) 


The basic equo^tions from (3.3a-e) for the plane case 
are as follows: 


5x 


V 

u 


(4.3a) 



(4.3b) 


0X 


+ 



(4.3c) 


u 



2 2'/ p ^ 

V + — f 1 



(im/ 

^Y+r 


2 

Y+1 


(4.4d) 


sin cr- 


(4.4e) 
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Sin cr and cos o~ are found in terms of 


s in 0 ” = 


Tji sf 


2 

cos cr 


f 


1/n 


+ sf 


where 


Br = 


2 „2/n , . 

n A and ii-j -j = 1 -n 


Boundar 3 r conditions follow from (2.4) 


P 


J2._ 

"f+l 


B, 


2 .,2 


+ sf 


V +1 
Y-f 


u = 1 •- 


2 

Y+1 


2- sf ijr^ 


Y+1 2/n Q o 

f + Bf r 


Bi t ijr 


1/n 


2/n p 9 

tf + sf 


y = 


follows; 


(4.5a) 

(4.51) 

(4.5c) 

(4.5d) 
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Out er E]cp aiisions 


X .and Y'j are the outer irdependent variahles' 

where 


X = X and (4.6a) 

w = ijr/x^ (4.61 d) 


'The transformed equations of (3.3a-e) v/ith the 
corresponding flov; •varie.hles replaced hy P, h, V, U and Y, 
are the following: 


8Y V ^ A w n 

^ = IT + -TTr-~' ^ 


(4.7a) 


§1 

0Vv' 


mr ^ 


(4.7h) 


IJ, - w n fer + 


9V 8P 

Qw 0W 


= 0 


(4.7c) 


2 2 2Y P 

U + Y + 


Y-1 D 


y-i/ 2 


= 1 


2(l-~) 


5/p = 


Br w 


Y+T i+l 


^2/n ^2(l-.n) 


(4.7d) 


C '• - ^2/n ^^i1-n)/n^2(^-n)" * -I 

(4.7e) 



42 


The shock conditions are transf omed into the following. 


+ 2/n 2(l-n)/n 2(l-n) 

A ' w X 


B, +. 

.4/n 2(l-n)/n 2(l-n) 

Aw X 


(4.8a) 

(4.8h) 


U = 1 


^2/n n)/n 


^^2(l-n)/n y2(l-n) 


+ » * . 1 


(4.8c) 


y+1 - [“J, _i) ^ ^ ■' 

A w X 




+. . .j 


(4.8d) 


Y = A w X 


(4.8e) 


evaluated at w = A. 


Outer Expansions 


The asymptotic expansions for the outer flow 


quantities can he directly written from the shock conditions, 


Pi(w) Pp 

P(X,w) = 

X^^1t 


(4.9a) 
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D(Z,w) 


D2(w) 

]3i (v/) + 

A 11 


4* 


(4.9b) 


Li (Zy V/ ) = 


U^(w) 

Z ^-1 1 


UgCw) 

^4nii 


•f 


(4.9c) 


V (X,v\/ ) 


V^(w) 

X^l l ?^11 


(4.9d) 


Y(X,w) = 


x”^ [ Y-] + 




• » • 


] 


(4.99) 


Substituting the above expansions into the equations 
(4.7a-e) we get the outer equations. Since the procedure is 
same as for the axi-symmetric case only the first order 
solution is given in this chapter. 


First Outer Equa,tion: 


= V., + 


A V/ n 



(4.10a) 


dY^ 

dv/“ 



— AV.jn.|.j ■ 


Awn 


dY. 


+ 


dP, 

W 


'■ = 0 


(4.10b) 


(4.10c) 


2U 


1 





(4.10d) 
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P. V 9 2(n-l)/n 

hr- = f 1-1 ) _Bl W 

'/+1 

(4.10e) 

jonditions 


P = 2„2 

(4.11a) 


(4.11h) 

rr 2 .2 2 

U = - A n 

(4.11c) 


(4.11d) 

Y = A 

(4. lie) 


Hure c^rain the outer solution is same as the self 
simili^r solution (plane case), the details of which are 
{.tivem in Appendix I. The first outer solution is the 
f ollov/ini,:. 

P,a,w) = aV hO) (4.- 12a) 




(4-12h) 


A^n [ fpx) + ^4-'^] - 


(-/+1 ) 


g(X) -• 

(4.12c) 


U^(X,w) = 1 - 
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a,«) = A n f(A) (4.124) 

(X,w) = aXx^ (4.12e) 

Holo.ci*jn botv'/een w and a follows from tlie eq_ual 7 ion 

(4.10a) 


[ (V+1)A- 2f(X) ] 


Imor Expansion: 

ijf and X are the independent variables and equations 
(^,3a-o) deocribo the inner flow field. Eor finding out the 
inner expansion, the same procedure as for the axi-symmetric 


adopted and the following expansions 

are arrived at. 

il 

2n^ A 

(4.13a) 

f(x,ijr) = 

?1 (t) 

(4.13b) 

u (x,f) = 

(f) 

^ -toVWTJZh ••• 

X 11 

(4.13c) 

V (x,f) = 

(t) , 

('/-23^^)/y 

(4.13d) 

7 = 

2n/Y 

(f) X ' + ... 

(4.13e) 
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Substituting these expansions in the equations (4.3a-e) we 
got the following inner equation 


2n^^y^ = '/v^ 

(4.14a) 

dT " “ 

. (4.14b) 

dp. 

af-= ° 

(4.14c) 

i Pi 

>^1 TV P1 = ° 

(4.14d) 

*’1 (4n/ 2 

Pit y+1 V+i T|-2/n_^ ^2 

(4.14e) 


The two Hatching conditions required to solve the inner 
equations are the following 

Pl (x,f ) = :^ n^ h ( 

Pron equation 4J4c 

P1 = 
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= S' (It) 


where 


I' (t) = h (\) H 


2B 


y^ = / F (f) df + A ??13 

^1 = ^ [ / !■ (t) dir + A ] 

_ -2 '■ A^ n^ h ( ^h ) 

^1 ( v+i )^ p ' (T|r) ' 

ConpoBiue Expansion: 

Pi fe, fr) = n2 h (X) -L_ 

I I 


o e. P = ^ g + [ !■ (T) - ^1 (f)] 


u 


1 p = 1 - tIt O') + 


Li_^ h (\) [ j l^r) “ iyirr ^ 


(V+1) 
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v. (x, If) = tV- a n f (>) [ / P (i|r) df - ifc /P(tJ)df ] 




(x, T^) = A Xx^ + 


[ / P (ljr)df- Lt / P (f) df ] 
l|F^ 


Irs^/i 


where P^ (f) = Lt P (ijl) 

The body for the case n = 2/5 is given in Pig. (14), The flow 
profiles are also given in Pigs. (15-1^ • 



CH APTER 5 

HIS CHS SICES AITD GOHC L HSIOIS 

Iho roGults of tlic first order axis3niiiiietric Hody for 
n in iho' r:mgo ‘2/3+o < n < 1 * are given in Figures (2-7). 
For the two cases j when n = 1/2 and n = 2/3, second order 
results arc given in Figures (2 - 3).’ Also the results 
of the first and second order solutions for the flow proper 
tics are given in Figures (8-13). In the case of the plane 
flow the results of the first order body and the flow 
properties are given in Figures (l4--16),for n = 2/3. Ratio 
of the body to the shock thickness for the axisyinnietric 
COG-.; is given in Figure (1?)* 

In Fifuirc (2), the first order results are the same 
ar. Ya’mra’s results for n = 1/2. Ihe difference, in the 
body shape by Sychov and Yakura, is due to the velocity- 
d-cfoct effect, vv.hich is included in Sychev's calculations 
but it; not included in the first order theory of Yakura. 
Th(; numerical relation obtained by Sychev is of a hi^er 
order th.an the first because he has used the two-term 
approximation for the velocity u and the one-term approxi- 
mation for the density. It is in close agreement with the 
present second order results. As n increases the body 
approaches the shock. Figure (17) shows that the entropy 
layer is of negligible importance when n > 0.65. In 
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Fifure (J) for n = 0,99> the body is almost strai^t making 
an ai'glc of 50.5^ A'hich is very nearly the same as the flow 
dcfloci;ion angls for a shock angle of 45°. This body coin- 
cides with Yakura’s bodj'' (producing a hyperboloidal shock) 
foT' largo distances. 

The first order pressure gradient is zero and the 
Second order pressure gradient exists only when the power- 
law exponent, n, is 1/2. In the outer region the first order 
and the second order densities coincide. Figures (10 and 11) 
show that the u velocity decreases as we go towards the body 
and tin. gradientrs are more predominent in the inner region. 
Figures (9 and 12) show the gradient of the entropy func- 
tion. As. wc go av'.'iy from the leading edge of the body the 
gradient occurs only near the body. The entropy layer thick- 
ness d'-ercases away from the leading edge but it exists 
,an i'oL' large distances. The second order results are not 
■v-;.;ry significant though the correct the body and the flow 
properties in the right direction. 

Tnc entropy layer and the Prandtl’s boundary layer 
have common features like the displacement thickness and 
zero transverse pressure gradients. The influence of the 
entropy layer varies with the body or shock geometry and 
is more predominent for n < n*. Hornung has discussed the 
range of n for the attached layer, the free layer and the 
blast wave region. The flow in the entropy layer has gone 
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1'- - -tr. shock and does not satisfy the hypersonic 
i;.Lena< b'.tiy iH.',|Uircments, ' As a consequence the entropy 
layo:r : '.ov, dote, not oodiihit similitude (6). But the self- 
oimil.-T' isolutionc become significant for larger distances 
i'v'orn fib., nose of the body, and also vvhen n is in the range, 
n* < n < 1. 




iijpg]pix I 

In f-encral, the solution of the partial differential 
equations governing a gas dynamic situation, can only he fotuid 
hy complicated numerical methods. But for certain gas motions, 
it turns out that those equations can, he reduced to ordinary 
differential equations and one can obtain exact solutions in 
closed form or approximate solutions throu^ the use of compa- 
ratively simple numerical integration methods. Such motions 
arv- call', d self-similar. 

The solution to the small disturbance equations of the 
steady tv/o dimensional or axisymmotric hypersonic flows come 
under i-hc class of similar solutions. The equations governing 
hypcrnonic flov; over slender two dimensional or axisymmetric 
body a:a. r-inilar to the unsteady hypersonic flow equations 
in one cp'ice dimension, the solution of which is known ( 4 ). 
This and the fact that the hypersonic motion of a slender body 
through gas produces only a transverse displacement of the 
gas (hypersonic equivalence principle) helps in assuming that 
the solution at one section (at a point on the body perpendicu- 
lar to the flow direction) is similar to the one at any other 


section on the axis. 
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Ono-dirensional IJnsteacly Equations: 

+ ly (pv) + 3y“^ = 0 


(A. 1a) 


^ 4. rr ^ JL ^ §£ 

ii?+^3y + r8y = 


0 


(A. It) 


(p/^/ ) + V ~ (p 4/) = 0 (A. 1c) 

j = 0 Plane flov</ 

j = 1 Axisymetric flow 

whore V)» r.nl v are the pressures^ density and velocity in 
tt‘, y direction rcspoctivoly. The equations A.1 can be 
trai 4 .sf i/fir.od, int tli<j hypersonic small disturbance equation 
by sii.n lt transf ornation x = U^t where is the free stream 
vcl-'city. In yonoral, Hayes transformation enables one to 
trannror..i the approximate equation representing a steady n- 
diricnsinnal flow into exact equation representing an unsteady 
n-l'-diricnsional flow. Thus all the results of the unsteady 
gr.sdynanics can be used to investigate the corresponding 
approximate steady flow. The time axis transforms into the 
X axis* This principle is called Hayes equivalence principle 


or the »Iaw of Plane Sections* 
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S lilt ion f the Ex.’^,ct Equations; 

jJi; ; ’nsionless variables V, Rj and P are defined to 
the ;'ep:.;ndcnt variables 


V == 

yt~’'v 

(A. 2a) 

r: 

Ay"^R 

(A. 2b) 

P = 

Ar" t ^P 

(A. 2c) 


whore 4 is a constant, the dinensions of which also contain 
the r.iasH since the dimensions of the quantities p and p also 
contain the mass; w is a dimensionless parameter. 

An additional parameter B is introduced so that its 
dii-n^-noions do not c 'ntain the mass, and number of independent 
viT'iabltc is na'uccd to one by defining a new independent 
varianleo ** , 

/ = Byt ^ ( 4 . 3 ) 

whore k A- dimensionless exponent. Ihe dimensionless variable 

Ic 

A is Oi-nEtant on a path for which y is proportional to t or 

k 

in an (.equivalent steady hypersonic flow to x . Substituting 
A, 2 and A. 3 in A, 1, we get the following set of ordinary 
diffcjcntial equations: 

>, [ (V-lc)7’ + J = -V(V-1) - (2-w)| 


(a, 4a) 
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.,[¥'+ (V-k)|l 1 = - (l+j_w)V (A.4b) 

[A(V-k) ] [ = -2(V-1) -w(Y-1)V 

(A. 4c) 


Hurc tho primes indicate differentiation with respect 
to .X . Body or a shock-wave in the flow field must lie on a 
liiiu of Oi nstant X , 

A now indcpt.ndent variable z is introduced to replace 
the v.ariahle P. 

2 = / (a. 5 ) 

ninrlv^ differential equation is obtained from A.4 
in whio’i P, R .■.nd oo not appear. 

[2 (v-i )+d ( '-1 )Vj (V-k)^+(Y-1 ) ( 1 -k)V (V-k)- [2 {V-1 )+E:('/-1 )z]] 

(V-k) I V(V-1)(7-k)-[(1+a)V-E]z] 

(a. 6) 



where K = [ 2(l-k) + kw ]// 

From the sane equation A. 5 following relations may be 


obtained • 


d log > 


T(V-1 ) (V-k)- [(1+3 )V-K]2 


(A.7a) 
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cl log(V-k )R _ __ (l+n--w )7 

d log A (7_l5-) . 


(A. 71)) 


Once a solution is otitained from l»6 for z(Y),X 
and R May bu obtained by quadratures over 7, 


Boundary Conditions: 

L sinilarity line is defined as the line in (y,t) 
space for which A is const'ant and y is proportional to t^. Hie 
velocity of a point on this lino is kyt"''. Brom A. 2a the velo- 
city of a fluid particle relative to a similarity line is 
(7-lc)yt , Thus 7=k is the boundary condition on a physical 
boundary which follows a similarity line. 

In the undisturbed region of flow the pressure must 
be consto.nt in both time and space. Iron A.2.c either k=2/(2-w) 
with P proportional to or P = 0. Except when k=2/(2-w) 

the pressure in the undisturbed region in front of a shock 
must be negligibly snail. Also density in the undisturbed 
region must be independent of t. Prom A. 2b it can be deduced 
th.'it R is a constant and the density is proportional to a power 
— w of y. In our problem the density is a constant in front of 
the shock and so w is taken to be zero. The shock conditioi^ 
give us the boundary values for the variables immediately 


behind the shock. 
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2k 

(A. 8a) 

Y+1 

(A, 8b) 

.y-1 

2k^ 

7?T 

(a. 8c) 

2/C/-1)k^ 

(Y+1) 2 

(A.8d) 


* In the solutionis of interest to us v\/hich contain the 
shockSjthe location of the hody is unknown in advance, The 
fact that this boundo-ry is a floating one requires that we 
inpose an additione,! condition to determine the solutions 


7- 


b 


k 


whore suffix b refers to the physical boimdary. 


Constant Energy Solution: 

Sodov"^ obtained unique integrals of equation A, 6 for 
the case of constant energy by considering conservation of 
energy across two similarity lines and X 2 * 

2 = (a-.1)7^(k-7)/2(7-k/Y) (A.9) 

Although the integral A, 9 is obtained from considera- 
tions entirely from those, we obtained equation A. 6, yet by 
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o.ctual substitution we find that A, 9 satisfies A, 6 exactly. 

For this solution k is given by 

k = 1/2+j-w (A. 10a) 

when w = 0, constant density case, 

k = 1/2+3 (A.IOb) 

The solution A. 9 satisfies the shock boundary condition 
A.8d. The quantity z is infinite at V=k/V . By actual computa- 
tion wG find that as A - 0, V - k/-/ and z - So we may set 
the inner boundary condition. 

\ = 0 5 V-. = k/y (A.il) 

b ’ b ^ 

This shov/s that at the core when a-j^= 0 or y=0 the density 
^ - 0 and the tenpero.ture T - while the pressure remains 
finite. The qurdratures of A. 7 together with the basic solu- 
tion A. 6 yield the results for a given value of time or x 
as the case may be 




7 


s 


= >. 



2k ^ 


f 


_ izl 

?s 


R = g 
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Ps 


V +1 

2k^ 


)2 


P = li 


I'lio functions given in ( 4 ) f, g, li are tabulated for 
the case of violent explosion when 3 = 0 or 1 in the following 
pages . 



TiiBIE 1 



O 

11 

il 

C\J 

ll 

A 

f 

g 

ll 

1 .0000 

1 .0000 

1.0000 

1 .0000 

0.9797 

0.9699 

0.8625 

0.9162 

0.9420 

0.9156 

0.6659 

0.7915 

0.9013 

0.8599 

0.5160 

0.6923 

0.8565 

0.8017 

0.3982 

0.6120 

0.8050 

0.7390 

0.3019 

0.5457 

0.7419 

0.6678 

0.2200 

0.4904 

0.7029 

0.6263 

0.1823 

0^4661 

0.6553 

0.5780 

0.1455 

0.4437 

0.5925 

0.5172 

0.1074 

0.4229 

0.5396 

0.4682 

0.0826 

0.4116 . 

0.4912 

0.4244 

0.0641 

0.4038 

0.4589 

0.3957 

0.0536 

0.4001 

0.4161 

0.3580 

0.0415 

0.3964 

0.3480 

0.2988 

0.0263 

0.3929 

0.2810 

0.2410 

0.0153 

0.3911 

0.2320 

0.1989 

0.0095 

0.3905 

0.1680 

. 0.1441 

0.0042 

0,3901 

0.1040 

0.0891 

0.0013 

0.3900 

0.0000 

0,0000 

0.0000 

0.3900 



TiiBEE 2 



j = 1 


f 

1 .0000 

1.0000 

0.9988 

0.9996 

0.9802 

0.9645 

0.9649 

0.9374 

0.9476 

0.9097 

0.9295 

0.8812 

0.9096 

0.8514 

0.8725 

0.7998 

0.8442 

0.7638 

0.8094 

0.7226 

0,7629 

0.6720 

C ».7242 

0.6327 

0 . 6894 

0.5989 

0.6390 

0.5521 

0.5745 

0,4943 

0.5180 

• 0.4448 

0.4748 

0.4073 

0.4222 

0.3621 

0.3654 

0.3133 

0.3000 

0.2571 

0.2500 

0.2143 

0.1500 

0.1286 

0,0000 

0.0000 


1.4 

k = 1/2 

g 

h . 

1 .0000 

1 .0000 

0.9973 

0.9985 

0.7653 

0.8659 

0.6285 

0.7832 

0.5164 

0.7124 

0.4234 

0.6514 

0.3451 

0.5983 

0.2427 

0.5266 

0.1892 

0.4884 

0.1414 

0.4545 

0.0975 

0.4242 

0.0718 

0.4074 

0.0595 

0.3969 

0.0362 

0.3867 

0.0208 

0,3794 

0.0123 

0.3760 

0.0079 

0.3746 

0.0044 

0.3737 

0.0021 

0.3733 

0.0008 

0.3730 

0.0003 

0.3729 

0.0000 

0.3729 

0.0000 

0.3729 
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: - SHOCK DISTANCE -IN ; AXiSYMMETRiC 
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